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EVALUATION 


The use of prior data in stru_ ‘uring reliability tests is attractive 
because of a potential for reducing test times (when the prior is 
favorable) and for structuring tests based on defined risks more 
meaningful than the classical producers and consumers risks. The 
technical basis for these tests is described in, RADC-TR-76-294, 
"Reliability Acceptance Sampling Plans Based Upon Prior Distribution". 
Unfortunately, determining the parameters of the prior and designing 
test plans based on these parameters have been formidable mathematical 
exercises. The object of this contract was to solve the latter 
difficulty by developing a simple method for designing tests from a 
given prior. Its product is a series of tables which will directly 
provide the test length and allowable number of failures from the 
parameters of the prior and the selected test risks. Although voluminous, 
the tables are easy to use and require no special mathematical skill or 
computer programming. Coupled with other studies recently completed to 
provide a simple means of determining the parameters of the prior, this 
report will be the basis for an RADC exhibit on reliability testing 
based upon prior distribution. 
ony Captor 


ANTHONY COPPOLA 
Project Engineer 
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1. INTRODUCTION 


This report presents a procedure for the design of reliability 
acceptance sampling plans for the case when the failure distribution 


is exponential, i.e. 


£(t|e) = 6 texp(-t/e) ; t>0, 86> 0 (1) 


where t denotes the time of failure and 6 is the mean time 
between failures (MTBF). The parameter @ is considered to be a 
random variable with known prior frequency distribution g(8). This 
prior distribution is assumed to be the conjugate inverted gamma 


density 


a. 


g(6) exp(-y/8)/T (A) ; Nadie 2 Ol (2) 


The successive lots or systems are taken to have different (con- 
stant but unknown) i=1,2,3,...; drawn independently from g(@). 
The procedure given here is applicable when accept/reject 

decisions are to be made regarding a sequence: vf systems or a 
sequence of lots. It is assumed that the decisions are based on 
data subject to random fluctuations, that a definite measure of 
loss is associated with each of the two decisions when they are 
inappropriately taken and that each lot or system is sentenced on 
its own merit without regard to the decisions taken about other 
lots or systems. 


Following four types of single sample plans are considered: 


(i) truncated plan with replacement, (ii) truncated plan without 





replacement, (iii) censored plan with replacement and (iv) censored 
plan without replacement. For a lot, which consists of a large 
number of components, all the four plans are feasible. For a system, 
which can be repaired as new upon failure, 100% inspection is assumed 
and plans without replacement are infeasible. For specified combina- 
tions of producer-consumer risks, the tables given here can be used 
for the design of a truncated single sample plan for a system. 

Other single sample plans can be obtained by using the equivalence 
relations given in this report. 

A large number of producer's and consumer's risk combinations 
can be obtained. However, in this report we provide tables for the 
relatively more useful combinations, viz (a,8*), (a*,8*) and (a,B). 
These risks are defined in Section 3. Tables for other risk com- 


binations can be obtained if needed. 
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2. EQUIVALENCE RELATIONS 


It is easy to note that plans with the same operating char- 





acteristic (0.C.) curve are equivalent i.e. plans with identical | 
P(A|6) vs @ curves are identical. We now give the conditions under ] 
which different types of single sample plans have the same operating | 
characteristic curve as the truncated plan for a system. That is, 

we give the equivalence relationships between a single sample 

truncated plan and the other types of plans. A detailed derivation 


of these relationships is given in Goel and Joglekar [3]. 
202 Plans With Replacement 


Zekek Truncated single sample plan for a system 


Under this scheme a single system is placed on a life test of 
duration T . Without loss of generality, repair is assumed to be 
instantaneous. If the observed number of failures, r, is less than 
or equal to a value r* i.e. r<r*, the system is accepted. Other- 
wise, it is rejected. Since the failure distribution is exponential, 
the number of failures in fixed time T has a Poisson distribution. 
Then, the probability of acceptance for a given 6 , can be written 


as: 


* -T/9 
P(alé) = r oS ry )® 


‘nt zi (201) 





alee Truncated plan for a lot 


In this scheme a sample of n items from a lot is placed on 
a life test of duration to - Failed items are immediately replaced. 
If r<r' , the lot is accepted, otherwise it is rejected. 


The failure process can be considered as a Poisson process 


with rate (n/é). Hence, 


“a -nt)/@ 


t : 
P(Aje,n) = P(x < r'|6.n) = £ i i sl (2.2) 


r=0 es 





It is easy to see that the 0.C. curves given by Equations (2.1) 

and (2.2) are identical if nt =T and r'‘=r*. Thus, the desired 
plan for this situation will be given by (t)=T/n, and r‘=r*). A 
Suitable value of n is chosen to make the sample representative 


of the lot and/or from cost considera ions. 


Zahs3 Censored plan for a system 


A system is life tested until r' failures occur. If the 
total test time t'>T', the system is accepted; otherwise it is 


rejected. If tyrtoreeertie represent the times to failure than 


yt 
Eg" = £ t; has a gamma distribution given by, 

i=l 

—t' 

= (r'-1) 

' 
£(t'|6) = £3+ t 
& re") 


Hence, 


-t' 

ae nr 1 

P (ale) = P(t" > T'I6) = f a tt) at! 
8 


T' 


tte") 











If t"=t'/o , then 


© ot" rial 5s ed | 2 NG 
Pile) = ¢ St d's t SUMAN (2.3) 
Tf ((r') r=0 r 


Equations (2.3) and (2.1) are identical if r'‘'=r*+l and T'=T. 


Hence the desired plan is (T,r*+l). 


2.1.4 Censored plan for a lot 


A sample of n items from a lot is life tested until r'‘<n 


failures occur. If the total test time t'>T' , the lot is accepted; 


otherwise it is rejected. 


' 
It is shown in Epstein (1960) that ae my . Hence, 
' ' 
p(ale) = P(t! > T' Je) = P(e > 22) 
8 0 
eo 1 = 
ws Le 1 W235, 
at’ P(r") 2 
F 
2t' 
where <a 
8 
If t"=>, then 
ba See oe r'-l Rad i 2 r 
Pale f Sit. gre t = (T'/8) docks 
F rte") r=0 r! 


Comparing Equations (2.1) and (2.4), the required plan is obtained by 


by taking r'‘'=r*+l and T'=T. 








Zee Plans Without Replacement 


waaak Truncated plan for a lot 


In this scheme a sample of n items from a lot is life tested 
for duration to Failed items are not replaced. If r<r' the 


lot is accepted; otherwise it is rejected. r has a binomial dis- 





tribution and 
j 

x" “t,/e ~t,/e | 

P(A|6@) = P(r ¢ r'[6) = £ () (l-e or (e of ner (2.5) 

r=0 | 

Epstein (1954) shows that for small a,8 and 8y/ty > 3 » @pprox- ; 
| 


imately r'=r* and 


G~ x 2r*) os | 
tot ot gn GCE). | 
2r* n 
Since 
2 
8X” 42, (2r*) 
UY et An 
2 
ft n 
we have ty) == an (~e) : 








eudee Censored plan for a lot 


n items from a lot are life tested until r'‘<n failures 
occur. Failed items are not replaced. If the total test time 
t'>T' , the lot is accepted; otherwise it is rejected. 

Epstein (1960) shows that, as in the replacement case, 
2t'/d v ae - Hence the designed plan is r'=r*+l and T'=T. 


The relationships between the above plans are summarized in 


Table 2.1. 





BQUIVALENCE RELATIONS 


Truncated Censored | 


Truncated Censored 
(T,r*) — (T,x* + 1)) (ty= T/n,r*) 
ze On GEog) rt) (Trt + 2) 


With Replacement (T,r* + 1 


Without Replacement 











3. RISK CRITERIA 


Since the parameter 6 is assumed to have a prior density 
g(a), the producer's and the consumer's risks may be quantified in 
various ways. The definitions of the various risks are now given 
in terms of P(A|8) and g(6). For the interrelationships and 
interpretations of various risks, the reader is referred to Goel 


and Joglekar [1976a). 


Sed Producer's Risks 


Classical: 


The probability of rejecting a system with 6 equal to the 


Specified value @) is equal to a 4 hats 
a = P(R|@), orl- a= P(A|@) (3.1) 
Average: 


The probability of rejecting a good system (6 > 84) is equal 


to a» 1.6. 


or 


(3.2) 


1-Ge - P(A|@) g(6)d9 
9. fF 9le)de 
% 








a en san amine es iieseinkss 4 . ace Sina anil 


tf 


Posterior: 
The probability of a rejected system being good is equal to 
a*, i.e. 


at = P(6>8.|R) | 


or 


g {1-P(A|8)} g(8)de 
at = 2 





= (3.3) 
1- sp (al 6)g(6)ae 
oO 


Probability of Rejection: 
The probability that a randomly selected system will be 


rejected is prespecified. This probability represents the long 
range fraction of the submitted systems that will be rejected by 
the sampling plan and is given by 


P(R) =1- f P(a}®) g(6)d8 (3.4) il 
0 i 


3.2 Consumer's Risks 


Classical: 
The probability of accepting a system with minimum acceptable 


MTBF 84 is equal to 8B, i.e. 





B = P(a|6,) (3.5) 


So 


Average: 


The probability of accepting a system of unacceptable reli-~ 


ability (8 < 9) is equal to 8» ist 


B= p(alée < 8) 








or 


(3.6) 


. 
f~* Pals) g(8)de 
0 


6 
Pt g(eyae 


Posterior: 
The probability of an accepted system being bad (6 < 81) is 


equal to s6* , i.e. 


Be = p(B < 6, la) 


or (3.7) 


a 
S°P(A| 8) g(6) de 
a) 





Be = 
Pp(a| 6) g(6)de 
Vv 


Alternate Posterior: 

To provide a better control on the distribution of 6 in 
the accepted systems, we define an additional posterior risk, g**, 
which specified the probability that an accepted system will have 


an MTBF less than the specified value 85 « Thus, 


Bet = pie < 8 lay 


or 
89 
f “P(al®) g(®)ae 
Bee = (3.8) 
f p(al8) g(8)ae 
0 
The various risk criteria have been defined above in terms of 
P(A|6) and g(6). Explicit expressions for the risks for all the 


four types of single sample plans can be obtained by substituting 








the corresponding expressions for P(A|6) derived earlier. 





Risk 


expressions for the truncated plan for a system are given below. 


ee Expressions for Producer's Risks 


Classical: 


e"/%o «rye yt 





xr 
l-a= J 
r=o Se! 
Average: 
ont oe ea. 
ge ee ee aU 
i = aa} =0 Yr. T (A) 
of —— gtAtl) .-¥/8 a6 
o [(A) 
Posterior: 
ts) r* -T,/'6 4 oN aa a 
prui-e i & (T/ 0} 7 gOrl, ¥/8a6 
1 « gue r=0 rs T (A) 
ie r* -T/6 a ae 
an 27 ¢ e (T/8) +x 8 (A+1) e~Y/ 84a 
° r=o Es Ta) 
Probability of Rejection: 
co =* bo i A 
Binke as Pps Sy Ee was 
° r=o0 x! T(A) 


3.4 Expressions for Consumer's Risks 


Classical: 


e tf S2 (T/01)* 


Ee 


(3.9) 


(3.10) 


(3.11) 


(3.22) 


(3.13) 





me 


—- 


Average: 
6 Ee F/G 2.8 <4 432) «6 
pos e 27) ee ae pe ¥/ as 
- BS oO Yr=o r! T(A) 
aS Se 
4 
° T (A) 
Posterior: 
6 © -T/0 r . 7 
pis SA per) Pig 


k= ¥=0 Ec T (A) 
Ss ee ene 


era) yy" pias) Wee a 
0 rs Tr (A) 


6 r x ne 
Bact oy gery ax78 a, 
gtt= OF rt P(A) 
el ee i: sc eee ne Gee aeerneanenees 
. 8 (r7o)* } * g{A+l) ev/? an 
° r=0 r! T(A) 


Oo 


12 





te Fig tape tht 


(3.14) 


(3.15) 


(3.16) 
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4. PDF AND CDF OF INVERTED GAMMA PRIOR 


Before designing the test plan, it is useful and instructive 
to study the probability density function (pdf) or the cumula- 
tive distribution function (cdf) of the inverted gamma prior distri- 
bution. For this purpose graphs and tables of this distribution 
have been provided in Appendix A for values of A from 0.4 to 4.0 in 
steps of 0.4 and for y*=y/8q from 0.4 to 4.0 in steps of 0.4. 

To see the usefulness of these graphs and tables, let us con- 


sider a few numerical examples. 


Example 1. Suppose from production data it has been determined that 
the prior distribution for 6 can be taken to be an inverted gamma. 
Let 85 be the specified value of @ and the discrimination ratio 
be equal to 2. Then 6* = 8/85 will also have an invertec gamma 


distribution. Let the pa._vaeters of this distribution be A = 2.0 


and y* = 1/85 = 0.4. From Table A.5-1 we find that 


P(e<0,) = P(o* < —l) = p(e* < 0.5) = 0.809 
and 
89 
P(0>6,) = P(o* > <2) ~ p(e* > 1.0) 1-P(8* < 1.0) 
0 


Values corresponding to these points are also shown in Figure A.5-l. 


1 ~ 0.938 = 0.062. 








i 


a 


- 


# 


What this information tells us is that a large percentage (80.9%) ot 
the area of the prior distribution lies below the minimum acceptable 
MTBF 8) while only 6.2% of the area lies above the specified MTBF 


8 In other words, the prior is not a good one and we need to 


0° 
use caution in accepting it as a prior before designing the test 


plan. 


- 


Example 2: For the same situation as in Example 1, let the para- 
meters be \ = 2.0 and y* = 2.0. Now, from Table A.5-5 and 


Figure A.5-5, we get 
PiGe -< 81) = 0.092 , PCS em) = .594 . 


This is the case of a good prior and we will need very little 


testing for acceptance purposes. 


Example 3: Now, let the prior parameters be A = 3.2, y* = 2.4. 


From Table A.8-6 and Figure A.8-6, we get 


P(@ < 84) = 0.169 and P(e > 85) = 0.382 


Here the prior seems reasonable and there is no reason for concern. 
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5. NUMERICAL PROCEDURE FOR THE DESIGN OF PLANS 


Plan design for a specified pair of producer-consumer risks 
is achieved by solving the corresponding pair of equations. Alter- 
nately, one may only consider the design of truncated plan for a 
system and obtain the remaining plans by using the equivalence 
relations given in Table 2.1. 

As an example, if the design criteria are P(R) and g* then 


given (¥r446y +6 q/P(R) -8*) Equations (3.12) and (3.15) can be solved 


to ‘ebtain Tf and r* . let T*=T/6 54 aaa K=6 5/8) and 
o=8/8> - Then Equations 3.12) and (3.15) become 
wo *® -T* Do » ik : 
rere te ne ee 38/0 reser _yy* 4-OFDET*/Fg9 (5.1) 
oO r=0 re T (A) 
and 
1/K r® -T*/¢ re A 3; 
meee Se CES AE Lg Vest ag 
ge= o.. eae = = - T (A) : (5.2) 
Wi { zr e (T*/o) ¥* ¢ (A+1)-y*, > 1 
0 r=9 ES P(X) 


Clearly, given (y*, %, K, P(R); @*) the values of © and r*® can 
be numerically obtained, if the solution exists. 

The numerical procedure consists of searching for those values 
of T and r* in the design region that satisfy the desired risk 
criteria. Let us call the consumer's risks g*, 8**, B as the left 
hand risks (LHR) and the producer's risks a*, a, P(R) as the right 


hand risks (RHR). It is intuitively clear that the LHR wil! increase 
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a as Fa a 


with increasing r* and decreasing T ; while the LHR will increase 
with decreasing r* and increasing T . Therefore, the region of 
feasible plans is a V-shaped area bounded by LHR und RHR as shown 
in Pigure 5.1. 

In order to search for the plan for the specified prior and 
specified risk criteria etc., the search region for T and r* 
is first delineated. We then set r* equal to its lower limit and 
use a binary search method to find a TL such that the plan 
(T»r*) yields a LHR closest to the desired value, the closeness 


being specified a-priori. For example, we can require 


LHR of (T..r*) - desired LHR 





<< e02 . 
Desired LHR 


For these values of (Ther*), the corr sponding RHR is computed. 

1f the computed value is less than or equal to the desired value 
then the plan found is the desired one. Otherwise we increase r* 
by 1 unit and continue the search by repeating the above procedure. 
The search is continued until we find a plan or we do not find one 
in the delineated region. In the latter case we increase the size 
of the region and proceed as discussed above. A similar procedure 
is followed if the RHR is first computed and the test for meeting 
the LHR criterion is satisfied next. The plans so obtained will 


not be exactly the same in general but will approximately satisfy 


the desired risks. 











\ RHR 
(gk, ar P(R)) 





‘Feasible Region 
for Plans 





of T of T 


Figure 5.1. Solution Regions For The Numerical 
Procedures For Plan Design 
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6. DESIGN OF TEST PLANS 


6.1 Design Tables 


Given (y*,\,K,a',8'), where y*=7/85 . K is the discrimina- 
tion ratio and a',8' are the desired producer's and consumer's 
risks, respectively, the problem is to obtain the design values T* 
and r* which satisfy the desired risks. Using the numerical 
procedure of Section 5 we have computed the values of T* and r* 
for several cases of interest for A from 0.05 to 4.0 in steps of 


0.05 and for \” from 0.0 to 3.45 in steps of 0.05. The various cases 


ares 

1. a= 6* = 0.10, K = 1.5 
2: a= §* = 6,16 , K = 2.0 
3. @ = §* = O10, K = 3.0 
4. a= B* = 0.15, K = 1.5 
5 a= p* = O15 ,; K = 2.0 
6. aa 6* = O.15 » K = 3.0 
7. o = §* = 6.20 | K 2 4.5 
8. a = 8* = 0.20, K = 2.0 
9. a = §* = 0.26, K = 3.0 
10. oa = 6* = 0.10, K= 1.5 
li. af = 6* = 0.10% K= 2.0 
12. a* = 6* = 0.10, K = 3.0 
13. at = p* = 0.15, K# 1.5 
14. at = 6* = 0.15 , K * 2.0 
15. o* @ 6% = 0.15, K = 3.0 
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16. a® = 6* = 6.20, K= 1.5 
17. a* = 6* = 0.20, K = 2.0 
18. a* = §* = 0.20, K = 3.0 
19. ©=8=0.10, K= 1.5 
20. a= 6 = 0.10, K = 2.0 
2. G8 Ge 0.258 . K = 1.5 
22. a@= 2 = 6.15 , K = 2.0 
23. ao = 6 = 0.20, K=1.5 
24. ao = B= 0.20, K = 2.0 


The design values for the above twenty-four cases are given 
in Tables 6.1 to 6.24 respectively. It should be pointed out that 
in obtaining these tables, the actual risks are less than or equal 
to the. specified values. In some cases the actual producer's risk 
will be close to the specified value while the consumer's risk will 
be much smaller than the specified value to get a feasible, minimum 
test time plan. 

Several other pertinent points about these tables are summarized 
below. 

(i) No plans exist for y* = 0. 

(ii) When a **** appears for the T* and/or r* values, the 


prior is not good and the test plan calls for rejection 


without testing. 
(iii) A plan T* = 0, r* = 0 means acceptance without testing. 
(iv) In cases where the prior is sharp (large \ values), there 
is a steeper change in plans than for the case of 


flatter priors (small A values). 


19 








(v) 


(vi) 


, 
} 
, 
| 
, 





In many cases the (T*-r*) surface, with respect to 


the risks, is quite flat. In such cases the computed 
test plans, due to numerical computations, can vary 
over a wide range If values without any significant 
effect on the producer's and the consumer's risks. 

In other words, different test plans will yield the 
same risks. 

The variation in the designed plans with respect to 

41 and y* values depends on the specific risk 
combinations. For (a, 8*) risks the plans follow 

a rising ridge from high y*, low A values to low y* | 
and high \) values as seen in Tables 6.1 to 6.9. For 

(a*, §®*) risks, the r* values first increase and then z | 
decrease as y* increases at constant A. When both | 
A and y* are varied simultaneously, the r* values ] 
first increase towards high \ and low y* and then | 
begin to decrease as seen in Tables 6.10 to 6.18. 

A similar pattern is seen for (a, 8)risks in Tables 


4 
{ 
6.19 to 6.24. | 
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6.2 Use of Design Tables 








We show the use of these tables via a few numericai examples. 





Example 1: 


Suppose the desired risks and the prior parameters are: 






= 


6.10, £46110, Ee 2 

















2.0 , ¥* = 0.4 








From Table 6.2 the plan is T* = 4.003, \’* = 3, i.e., conduct 


the test for 4.003 (89) time units and if the number of failures 





does not exceed 3, accept the system. 







Example 2: 





£0.10, #* «0.10, K= 2 


e| 
A 







2.0 






, 









From Table 6.2, T* = 0, r* = 0, i.e. the prior is very good and we 





can accept the system without testing. 







Example 3: 








From Table 6.11, the plan is T* = 1.281, r* = 2. 
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7. CONCLUSIONS 


In this report we have given tables for the design of fixed 





time reliability test plans based upon a prior distribution for 6 





A total of 24 sets of tables were included for various values of 





(a,8*), (a*,8*), (a,8) and discrimation ratios 1.5, 2.0 and 3.0. | 
| 
Also given are plots and tables for the pdf and cdf of the inverted 


gamma distribution, the conjugate prior for 6 


| 
. 
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APPENDIX A 
TABLES AND GRAPHS OF PDF ANU CDF OF 


INVERTED GAMMA DISTRIBUTION 


In this Appendix tables and graphs are given for the pdf 
and cdf of inverted gamma distribution for A from 0.4 to 4.0 
in steps of 0.4 and for y* = ¥/85 from 0.4 to 4.0 in steps of 


0.4. 





TABLE A.1-1 


Cumulati.2 Distribution Function of inverted 


Gamma Distribution with A = .4, y* = ¥/8q = .4 
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FIGURE 12 Ptots oF THe PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 = +4, 7¥ = <8. 
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TABLE A.1-3 


Cumulative Distribution Function of Inverted f 
Gamma Distribution with 4 =.4, y* = ¥/89 = 1.2 
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FIGURE *°  Ptots oF THe PDF anp CDF oF THE INVERTED | 
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TABLE A.1-4 


Cumulative Distribution Function of Inverted 
Gamma Distribution with »’4 =.4, y* = ¥/8y = 1.6 
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FIGURE 1-4 — PLots oF THE PDF AnD CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 = .4 , y = 1.6, 
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Cumulative Distribution Function of Inverted 
Gamma Distribution with }4 =.4, y* = 1/85 = 2).0 
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FIGURE 1-5 Prots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH } = -4 . Y = 2.0, 
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Cumulative Distribution Function of Inverted 
Gamma Distribution with A =.4, y* = ¥/8p = 2.4 
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FIGURE 1-6 PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH \ = .4, y =2.4 
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TABLE A.1-7 


Gamma Distribution with A’ =.4, 
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FIGURE 1-7 PLots OF THE PDF ano CDF oF THE INYERTED 
GAMMA DISTRIBUTION WITH 4 = .4 , Y =2.8 , 
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TABLE A.1-8 


Cumulative Distribution Function of Inverted 


Gamma Distribution with A =.4, y* = ¥/8 =e 
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FIGURE 1-8  PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =.4 , Y = 3.2, 
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Cumulative Distribution Function of Inverted 


Gamma Distribution with \} =.4, y* = ¥/85 = 3.6 
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TABLE A.1-10 


Cumulative Distribution Function of Inverted 


Gamma Distribution with A =.4, y* = ¥/85 = 4.0 
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FIGURE 1-10 PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =.4 , y =4.0, 
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TABLE A.2-1 


Cumulative Distribution Function of Inverted 
Gamma Distribution with 4 =.8, y* = ¥/85 = .4 
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FIGURE 2-1 PLots oF THe PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 = .8, ¥Y =.4 , 
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TABLE A.2-2 


Cumulative Distribution Function of Inverted 
Gamma Distribution with A =.8, y* = ¥/85 ci 
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TABLE A.2.3 


Cumulative Distribution Function of Inverted 


Gamma Distribution with } =.8, y* = ¥/8 9 = 1.2 
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TABLE A.2-4 


Cumulative Distribution Function of Inverted 
Gamma Distribution with X} =.8, y* = ¥/86 = 1.6 












0.000 3,500 0.527 8.500 0.740 











0.565 





0.005 





4.000 9.000 









0.500 0.026 4.500 0.597 9. 500 0.760 





0.624 





10.000 





0.062 5.000 








0.647 11.000 





0.082 53500 








0.145 6.000 0.668 12.000 












0.686 








0.261 6.500 13.000 











2.000 0.353 7.000 0.701 14.000 0.820 








0.425 7.500 0.716 15.000 0.839 












0.481 8.000 0.728 16.000 0.837 





MODE = .8889 


A-27 





SS 
a 


-00 


24,00 





OD 

i = we Lisi. 
«Oo a 

(om - 

Oo. Ry: 
6 e 
° ‘ee | 
rel a 

F 7 

= = 

«xe an 

a See 
0. 6.00 9.00 Le2vo 15.00 

THETASTAR 


FIGURE 2-4 PLoTs OF THE PDF anp CDF oF THE INYERTED 
GAMMA DISTRIBUTION WITH 1 =.8 , y =1.6, 


A-28 








- - Soar sere ~ — 
ae EE Ee 


TABLE A.2-5 


Cumulative Distribution Function of Inverted 
Gamma Distribution with \} =.8, y* = ¥/85 = 2.0 
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FIGURE 2-5 PLoTS OF THE PDF ANd CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 1 =.8 , y =2.0, 
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TABLE A. 2-6 


Cumulative Distribution Function of Inverted 
Gamma Distribution with } =.8, y* = y/8) = 2.4 
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Plots oF THE PDF AND CDF OF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 = -8, Y =2-4, 
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Cumulative Distribution Function of Inverted 


Gamma Distribution with A’ =.8, y* = ¥/85 = 259 
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FIGURE 2-7 PLots oF THE PDF ann CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =.g , Y =2.8 4 
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TABLE A.2-8 


Cumulative Distribution Function of Inverted 
Gamma Distribution with } =.8, y* = ¥/85 = 3.2 
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Cumulative Distribution Function of Inverted 


Gamma Distribution with 4} =.8, y* = ¥/85 = 3.6 
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FIGURE 2-9 PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 1 =.8 + Y =3.6 4 
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TABLE A. 2-10 





Cumulative Distribution Function of Inverted 
Gamma Distribution with \A =.8, y* = ¥/8o = 4.0 
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PLots of THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =-8 , Y = 4.0. 
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TABLE A.3-1 


Cumulative Distribution Function of Inverted 
Gamma Distribution with \ =1.2, y* = ¥/89 = .4 
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FIGURE 3-1  PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH » =1.2 4 Y = .4 4 
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TABLE A. 3~2 


Cumulative Distribution Function of Inverted 


Gamma Distribution with A =1.2, y* = ¥/89 = .8 
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PLots OF THE PDF AnD CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =1.2 - Y = .8 
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TABLE A. 3-3 


Cumulative Distribution Function of Inverted 
Gamma Distribution with A =1.2, y* = ¥/89 = 1.2 
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FIGURE 3-3  PLoTs oF THE PDF anp CDF OF THE INVERTED 
GAMMA DISTRIBUTION WITH % =1.2 4 Y *1.2 . 
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Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =1.2, y* = ¥/85 = 1.6 
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FIGURE 3-4 PLoTs OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH > =1.2 4 Y = 1.6 
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TABLE A. 3-5 


Cumulative Distribution Function of Inverted 


Gamma Distribution with A =1.2, -7* = ¥/8o NO) 
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FIGURE 35 PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH } =1.2 , y =2.0 
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TABLE A. 3-6 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =1.2, y* = ¥/89 = 2.4 
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FIGURE 3-6 PLots oF THe PDF AND CDF oF THE INYERTED 
GAMMA DISTRIBUTION WITH % =1.2 , Y = 2.4. 
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Gamma Distribution with A =1.2, y* 
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Plots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =1.2, ¥ = 2.8, 
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TABLE A. 3-8 


Cumulative Distribution Function of Inverted 
Gamma Distribution with \} =1.2, y* = ¥/85 = 3.2 
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FIGURE 3-8  PLoTs oF THE PDF AND CDF oF THE INVERTED 


GAMMA DISTRIBUTION WITH * =1.2 4 Y = 3.2. 
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Cumulative Distribution Function of Inverted 


Gamma Distribution with A =1.2 y* = ¥/8 9 = 3.6 
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FIGURE 39 Pots oF THE PDF Anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH } =1.2 , Y = 3.6. 
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PLots OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH * = 1.24 Y = 4.05 
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TABLE A.4-1 


Cumulative Distribution Function of Inverted 
Gamma Distribution with i =1.6, y* = ¥/8 = .4 
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FIGURE 4-1 | PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH * = 1.64 Y =.4 6 
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TABLE A.4-2 


Cumulative Distribution Function of Inverted 
Gamma Distribution with 4 =1.6, y* = y/8, = -8 
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FIGURE 42 = PLoTs oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH X =1.6, Y = .8, 
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Cumulative Distribution Function of Inverted 
Gamma Distribution with \ =1.6, y* = ¥/® = £2 
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FIGURE 4-3 PLots oF THE PDF anp CDF oF THE INYERTED 
GAMMA DISTRIBUTION WITH 4 =1.6, y = 1.2. 
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Gamma Distribution with A =1.6, y* 
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FIGURE 4-4 PLOTS oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =1.6 , Y = 1.6. 
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Cumulative Distribution Function of Inverted 


Gamma Distribution with A =1.6, y* 
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FIGURE 5-4 PLots OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH } =2.0, y =1.6 
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FIGURE 7-5 Ptots oF THE PDF AnD CDF OF THE INVERTED 
GAMMA DISTRIBUTION WITH } =2-8, y = 2-0, 
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TABLE A.7-6 


Gamma Distribution with \ =2.8, 


VARIANCE = 2.222 
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FIGURE 7-6  PLoTs oF THE PDF AND CDF oF THE INYERTED 
GAMMA DISTRIBUTION WITH * =2-8, y =?-4, 
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TABLE A.7-7 


Cumulative Distribution Function of Inverted 


Gamma Distribution with A =2.8, y* 


MEAN = 1.556 
VARIANCE = 3.025 


MODE = .7368 
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FIGURE 7-7 PLoTs oF THE PDF anp CDF oF THE INVERT 
GAMMA DISTRIBUTION WITH 4 = 2.8, 
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TABLE A.7-8 


Cumulative Distribution Function of Inverted 
Gamma Distribution with i =2.8, y* = ¥/8 9 = 3.2 
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MEAN = 1.778 
VARIANCE = 3.951 
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PLoTsS OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH \ =2.8, yY =3.2 
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TABLE A.7-9 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =2.8, y* = ¥/89 = 3.6 
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0.917 9.000 
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VARIANCE = 5.0 


MODE = .9474 
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FIGURE 7-9 = PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 = 2-8, y = 3.6, 
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TABLE A.7-10 


Cumulative Distribution Function of Inverted 
Gamma Distribution with \ =2.8, y* = ¥/85 = 4.0 


0.863 8.500 
0.896 9.000 
0.919 9.500 
0.936 10.000 
0.948 11.000 
0.958 12.000 
0.965 13.000 
0.971 14.000 
0.975 15.000 
0.979 16.000 
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VARIANCE = 6.173 


MODE = 1.053 
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TABLE A.8-1 


Cumulative Distribution Function of Inverted 


MEAN = .1818 
VARIANCE = 2.755 


MODE = 9.524 





Gamma Distribution with i =3.2, y* 
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FIGURE 8-1 PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =3.2, Y =.4 , 
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TABLE A.8-2 


Cumulative Distribution Function of Inverted 
Gamma Distribution with A =3.2, y* = ¥/85 = .8 


0.999 8.500 
0.999 9.000 
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1.000 10.000 


1.000 11.000 
1.000 12.000 
1.000 13.000 
1.000 14.000 
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VARIANCE = .1102 


MODE = .1905 
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FIGURE 8-2 PLots oF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =3-2, y = -8, 
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TABLE A.8-3 


Cumulative Distribution Function of Inverted 
Gamma Distribution with } =3.2, y* = y/@, = 1.2 
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FIGURE 8-3 PLoTs OF THE PDF AnD CDF oF THE INVERTED 
* 
GAMMA DISTRIBUTION WITH 4 =3.2, Y =1.2 


A-146 





Aad eS eee Rees 


TABLE A.8~4 


Cumulative Distribution Function of Inverted 


Gamma- Distribution with \ =3.2, y* = ¥/85 = 1.6 





02993 8.500 

0.995 9.000 
0.500 0.426 4.500 0.996 9.500 1.000 | 
0.667 0.618 5.000 0.997 10.000 1.000 
0.750 0.687 5.500 0.998 11.000 1.000 
1.000 0.820 6.000 0.998 12.000 1.000 
1.500 0.928 6.500 0.999 13.000 1.000 
2.000 0.965 7.000 0.999 14.000 1.000 
2.500 0.981 7.500 0.999 15.000 1.000 
3.000 0.988 8.000 0.999 16.000 1.000 





MEAN = .7273 


VARIANCE = .4408 


MODE = .3810 
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FIGURE 8-4  PLoTs oF THE PDF anp CDF of THE INVERTED 
GAMMA DISTRIBUTION WITH \ =3.2, Y =1.6, 
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TABLE A.8-5 





Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.2, y* 
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FIGURE 8-5  PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =3.2 , yY =2.0 
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TABLE A.8-6 


Cumulative Distribution Function of Inverted 
Gamma Distribution with i =3.2, y* = ¥/9p = 2.4 


0.977 8.500 . 


0.984 9.000 
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0.993 11.000 
0.995 12.000 
0.996 13.000 
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MEAN = 1.091 
VARIANCE = .9917 
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FIGURE 8-6 PLoTs oF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =3.2, ¥ = 2.4. 
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TABLE A.8-7 


Cumulative Distribution Function of Inverted 


Gamma Distribution with A =3.2, y* = ¥/8q = 2.3 







8.500 
0.013 4.000 0.976 9.000 0.998 
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0.667 0.245 5.000 0.987 10.000 0.998 
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1.000 0.518 6.000 0.992 12.000 0.999 
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FIGURE 8-7 = PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =3.2, Y =2.8, 
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TABLE 3-8 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.2, y* = ¥/85 = 3.2 


0.000 0.951 8.500 


0.005 0.965 9.000 
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0.820 0.993 14.000 
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FIGURE 8-8 PLots OF THE PDF and CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH A =3.2 4 Y =3.2 4 
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TABLE A.8-9 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.2, y* = ¥/8q = 3.6 
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FIGURE 8-9 PLots OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH A = 3.2, Y =3.6, 
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TABLE A.8-10 


Cumulative Distribution Function of Inverted 
Gamma Distribution with \ =3.2, y* = ¥/8q = 4.0 
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FIGURE 8-10 PLoTs oF THE PDF ann CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH A =3.2, Y = 4.0, 
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TABLE A.9-1 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.6, y* = ¥/8o = .4 
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FIGURE 9-1  PLots oF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =3.6 4 Y =-4 + 
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TABLE A. 9-2 


Cumulative Distribution Function of Inverted 
¥/85 = .8 


Gamma Distribution with A =3.6 y* = 
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FIGURE 9-2 PLOTS OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH } =3.6, 7 = .8 
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TABLE A.9-3 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.6, y* = ¥/8 = 1.2 
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A-165 





De. Actas Bi deci entrain et ate Raa Tae 


FIGURE 9-3 PLoTs OF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4.=3.6 », Y = 1.2, 
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Cumulative Distribution Function of Inverted 
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Gamma Distribution with A =3.6, y* 
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FIGURE 9-4 PLoTS OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH  =3.6 + Y =1.6.4 
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TABLE A.9-5 


Cumulative Distribution Function of Inverted 
Gamma Distribution with A =3.6, y* = ¥/8y = 2.0 
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FIGURE 9-5  PLoTs oF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH } =3.6 - Y =2.0. 
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TABLE A. 9-6 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.6, y* 
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FIGURE 9-6  PLots oF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH \ =3.6. Y =2.4,. 
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TABLE A.9-7 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.6, y* = 1/85 = 2.93 
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PLOTS OF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH \ =3.6 , Y =2.8, 
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TABLE A.9-8 


Cumulative Distribution Function of Inverted 


Gamma Distribution with i =3.6, y* = ¥/8 9 = 3.2 
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FIGURE 9-8 PLots oF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH \} =3.6, Y = 3.2, 
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TABLE A. 9-9 


Cumulative Distribution Function of Inverted 


Gamma Distribution with i =3.6, y* = ¥/89 = 3.6 
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FIGURE 9-9 PLots oF THE PDF AND CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH \ =3.6 ». Y = 3.6. 


A-178 








TABLE A.9-10 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =3.6, y* = ¥/85 = 4.0 


0.200 0.000 0.950 8.500 


0.333 0.001 0.965 9.000 
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FIGURE 9-10 PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =3.6 + Y 4.0 . 
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TABLE A.10-1 


Cumulative Distribution Function of Inverted 
= ¥/8y = .4 


Gamma Distribution with \ =4.0, y* 
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TABLE A.10~2 


Cumulative Distribution Function of Inverted 
Gamma Distribution with A =4.0, y* = 1/89 = .8 


0.433 1.000 8.500 
0.778 1.000 9.000 
0.921 1.000 9.500 
0.966 1.000 10.000 
0.977 1.000 11.000 
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0.998 1.000 13.000 
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TABLE A.10-3 


Cumulative Distribution Function of Inverted 


Gamma Distribution with \ =4.0, y* = ¥/8 ere 
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FIGURE 10-3 PLoTs oF THE PDF AND CDF OF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =4.0 , Y = 1.26 
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FIGURE 10-7 PLots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH = 4.0, y =2.8, 
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FIGURE 10-8 Plots oF THE PDF anp CDF oF THE INVERTED 
GAMMA DISTRIBUTION WITH 4 =4.0, y = 3.2, 
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FIGURE 10-9 PLots oF THE PDF AND CDF oF THE INVERTED 
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GAMMA DISTRIBUTION WITH 4 =4.0, Y =4.0. 
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